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ABSTRACT 

We present the discovery of the first T dwarf + white dwarf binary system 
LSPM 1459+0857 AB, confirmed through common proper motion and spectroscopy. 
The white dwarf is a high proper motion object from the LSPM catalogue that we 
confirm spectroscopically to be a relatively cool (T ft = 5535±45K) and magnetic 
(B^2MG) hydrogen-rich white dwarf, with an age of at least 4.8 Gyrs. The T dwarf 
is a recent discovery from the UKIRT Infrared Deep Sky Survey (ULAS 1459+0857), 
and has a spectral type of T4.5+0.5 and a distance in the range 43-69pc. With an 
age constraint (inferred from the white dwarf) of >4.8Gyrs we estimate T e g = 1200- 
1500K and logg =5.4-5.5 for ULAS 1459+0857, making it a benchmark T dwarf with 
well constrained surface gravity. We also compare the T dwarf spectra with the latest 
LYON group atmospheric model predictions, which despite some shortcomings are in 
general agreement with the observed properties of ULAS 1459+0857. The separation 
of the binary components (16,500-26,500 AU, or 365 arcseconds on the sky) is con- 
sistent with an evolved version of the more common brown dwarf + main-sequence 
binary systems now known, and although the system has a wide separation, it is shown 
to be statistically robust as a non spurious association. The observed colours of the 
T dwarf show that it is relatively bright in the z— band compared to other T dwarfs 
of similar type, and further investigation is warranted to explore the possibility that 
this could be a more generic indicator of older T dwarfs. Future observations of this 
binary system will provide even stronger constraints on the T dwarf properties, and 
additional systems will combine to give a more comprehensively robust test of the 
model atmospheres in this temperature regime. 
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1 INTRODUCTION 

Large scale Near Infrared (NIR) and optical surveys such 
as the 2-Micron All Sky Survey (2MASS), the Sloan Dig- 
ital Sky Survey (SDSS), and the UKIRT Infrared Deep 
Sky Survey (UKIDSS) are aiding the identification of a 
rapidly increasing number of 'field' brown dwarfs (BDs) (e.g. 
Lodieu et all |2007| . iPinfield et all |200S| , iBurningham et all 
20091 ). as well as probing down i nto new cooler tempera- 
ture regimes (|Warren et all 120071 . IBurningham et al 1 l200Sl . 



* E-mail:adjoncs@das.uchile.cl. Based on observations made with 
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Delorme et ai]|2008l . IBurningham et al.ll2009l , iLeggett et ail 
2009h . In general the estimation of properties of these BDs 



(e.g. age, mass, metallicity) currently relies on model fit- 
ting. However, the models are very sensitive to a variety 
of poorly understood processes in B P atmospheres, su ch as 
the formation of dust conde nsates jAllard et al l200lh and 
non-equilibrium chemistry l|Saumon et al. I l2007h , and the 
spectroscopic fitting of atmospheric properties (T e g, logg, 
[M/H]) is a major challenge. Crucially, the nature of BD 
evolution means that the mass-luminosity relation depends 
strongly on age, and in the absence of well constrained atmo- 
spheric properties there is no way to accurately determine 
mass and age. 
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Identifying objects where one can pin down these prop- 
erties independently can help aid the calibration of mod- 
els. BDs that are members of binaries, where the primary 
member h as age constraints are good sources of benchmark 



member nas age constraints are good sources ot benchmark 
BDs (e.g. iDav- Jones et al.1 l2008j IBurningham et all [20091 . 



IZhang et alJl20ldlFahertv et alJl2009l ). In particular white 
dwarf primaries can provide a hard lower limit on the age 
of the system (from the white dwarf cooling age) and in 
the case of high mass white dwarfs (where the main se- 
quence progenitor star will have a short lifetime and the 
age will thus be essentially the same as the cooling age of 
the w hite dwarf) , could p rovide ages constrained at the 10% 
level (jPinfield et alj |2006 and reference therein). 

Binary systems containing a white dwarf and a 
BD however are observationally rare, and only a 
handful of such binaries have been identified. To 
date there are only five L dwarf + white dwarf 
systems that h ave been spectroscopically confirmed, 
GD 16 5B (L4; IZuckerman fc Beckljnl Il992l). GD 140 
(L6/7; iFarihi fc Christopher! |2004 Dobbie et alj 120051) 
WD 137 - 349 (L8; iMaxted et all 20061. [Burleigh et al.1 
20061) . PG1234+482 (LP; ISteele et al.1 l200l iMullallv et al.1 



20071 ) and PHL 5038B fL8; ISteele et alj|2009h previously the 



latest type BD companion to a white dwarf. There have also 
been several BD companions to whit e dwarfs found as cat - 
aclysmic variable systems (CVs, e.g. iLittlefair et al.l 120081 ) 
, although the nature of such objects means that they my 
be less useful in the context of studies of typical BD atmo- 
spheres. 

We present here results from our ongoing search to iden- 
tify widely separated BD companions to white dwarfs, ex- 
pandin g on our earlier search es of 2MASS and SuperCOS- 
MOS ( Day- Jones et al.l |2008f) , to include the first results 
from our combined search of UKIDSS and SuperCOSMOS. 
We present here the discovery of the first T dwarf + white 
dwarf binary system, which we confirm through common 
proper motion and spectroscopy. 

In Section [2] we describe the ongoing search to identify 
and spectroscopically confirm late T dwarfs in the UKIDSS 
Large Area Survey (LAS). In section [3] we describe our 
proper motion measurements of these confirmed T dwarfs 
and our techniques to search for binary companions to these 
objects. In section|3]we describe the spectroscopy of a white 
dwarf candidate companion and the resulting constraints on 
its properties. In section [5] we statistically assess the likeli- 
hood that our new T dwarf + white dwarf binary system 
is spurious. Section [6] determines constraints for the atmo- 
spheric properties of the T dwarf, taking advantage of its 
benchmark age constraints (from the white dwarf primary) . 
We also perform some basic spectral model fits to the T 
dwarf spectrum, and compare the resulting predictions. Fi- 
nally in Section[7]we conclude with further discussion of the 
system as a useful benchmark and comment on the direction 
of future work. 



2 FINDING T DWARFS IN THE UKIDSS 
LARGE AREA SURVEY 

The UKIDSS LAS has been searched for T dwarfs using se- 
lection techniques based on the observed UKIDSS+SDSS 
colours of previously identified T dwarfs, as well as the- 



oretical predictions for the cooler T e g= 400-700K regime. 
The selection techniqu es used to iden t ify th ese T dwarfs are 
described in detail in iPinfield et ail (2008). In this paper 
we consider T dwarfs that were spectroscopically confirmed 
(Burningham priva te comm.) by the Su mmer of 2008 (see 



Pinfield et al.1120081 and a sub-sample of IBurningham et al.l 
20101) . and the sky coverage appropriate to this sample in- 



cludes the full LAS second data release, 72% of the new sky 
in the third data release, and 66% of the new sky in the 
fourth data release. In total, this sample was identified in 
890 square degrees of LAS sky. 



2.1 Follow-up photometry and spectroscopy 

Candidates are followed up in general with additional imag- 
ing in the near infrared and/or optical. This allows con- 
firmation of the expected T-like colours and rules out var- 
ious forms of contamination (e.g. faint M dwarfs with low 
signal-to-noise and blue- ward scattered near-infrared colours 
due to large photometric uncertainty, as well as solar sys- 
tem objects that can appear as non detections due to their 
motion). This followup has been performed using a variety 
of telescope/instruments, including the Wide Field Cam- 
era (WFCAM) and Fast Track Imager (UFTI) on UKIRT, 
and the Long-slit Intermediate Resolution Infrared Spectro- 
graph (LIRIS) on the William Herschel Telescope (WHT) 
(all near infrared), as well as the ESO Multi-Mode Instru- 
ment (EMMI) and the ESO Faint Object Spectrograph and 
Camera (EFOSC2) on the New Technology Telescope (opti- 
cal). Spectroscopic confirmation of the LAS candidates was 
also achieved using a number of facilities, including the Near 
Infrared Camera and Spectrograph (NIRI) and the Gemini 
Near Infrared Spectrograph (GNIRS) on the Gemini tele- 
scopes, and the Infrared Camera Spectrograph (IRCS) on 
Subaru. The Near Infrared Camera Spectrograph (NICS) 
on the Telescopio Nazionale Galileo (TNG) and the UKIRT 
imager Spectrograph (UIST) were also used for the brighter 
T dwarfs. Details of the follow-up imaging and spectroscopic 
strategies, as well as the relevant redu ction and calibratio n 
techniques used are f urthe r discussed in lLodieu et ail [[20071 ) . 
IPinfield et ail (|2008l ) and IBurningham et al.1 i|2010l) . Corre- 
sponding spectral types were derived using the uni fied T 
dwarf classificatio n scheme of Burgasser et al.l ([2006 ) with 
an extension from IBurningham et al.l (|2008l ) for the latest 
types. 



2.2 T dwarf distances 

T dwarf distances have been estimate d using the abs olute 
magnitude-spectral type relations from lLiu et al.l (120061 ). as- 
suming the T dwarfs are single objects. We calculated Mj, 
choosing the J— band magnitude over the H— and K— 
bands, as models have suggested that the J— band may be 
less sensitive to variations in metallicity and gr avity than 
the H— and K— bands (e.g. iMarlev et all [2003 ). The un- 
certainties in the distance were obtained by taking into ac- 
count the error in the spectral type (ty pically ±0 - 5) an d 
the residuals of the polynomial fits from I Liu et al.l (2006). 
Fig. [T] shows the spectral type distance distribution for the 
spectroscopically confirmed sample of 49 T dwarfs that we 
consider in this work. It can be seen that the sample spans 



Discovery of a T dwarf + white dwarf binary 3 



100 




I , , , i , , , i , , , i , , , i , , , I 

2 4 6 8 10 

T dwarf spectral type 

Figure 1. T dwarf distance vs spectral type for our spectroscop- 
ically confirmed T dwarfs with second epoch imaging. 

the spectral type and distance range T2-9 and 12-80pc. This 
is comprised of two T2-3 dwarfs, 29 in the T4-6 range, six- 
teen T6-8 dwarfs and two T8+ dwarfs. Due to the cooler 
temperatures and thus fainter nature of later T dwarfs it 
can be seen that we are more sensitive to earlier T dwarfs 
out to further distances. 



3 T DWARF PROPER MOTIONS 

The photometric follow-up program provided second epoch 
imaging for the LAS T dwarfs, which we combined with 
the LAS images to give two epochs to calculate proper mo- 
tion. We used the iraf routines GEOMAP to derive a geo- 
metric transformation between the two epoch images, and 
GEOXYTRAN to apply these transforms to the T dwarf po- 
sitions. Centroiding uncertainties were calculated based on 
simulated data with appropriate Poisson noise injected. The 
availability and quality of the measured proper motions of 
the T dwarfs is dependent on several factors including: the 
baseline between the epochs, the number of stars that can be 
used for positional reference in each of the images, the S/N 
of both the T dwarf and the reference stars, and the prox- 
imity of the T dwarfs to the edge of the WFCAM detector 
array (in the first epoch images). We were able to measure 
the proper motions of 19 T dwarfs from the 49 strong sample 
of spectroscopic confirmations that we consider, using an av- 
erage of 12 reference stars across a baseline of 0.5-1.5 years. 
A vector-point diagram of the T dwarf proper motions is 
shown in Fig [21 where two T dwarfs found to have common 
proper motion companions (see Section ^. ip are highlighted. 

3.1 A search for common proper motion 
companions 

We searched for candidate common proper motion com- 
panions to our sample of T dwarfs with reliable proper 
motions, by selecting a magnitude limited sample of Su- 
perCOSMOS sources, where R<21 (the magnitude limit 
at which proper motions are measured) and have accu- 
rately measured proper motions such that PM/apM J?3. 



1 i i 1 1 1 i i 


1 

ULAS J 1459 + 0857 ~ 
* 


T 

: + * 

I Wolf 940B H ~^~ H 
i , , , 


_4Jt l ~^ H 

+ ; 



1000 500 -500 

Proper mo;ion RA (mos/yr) 

Figure 2. A vector point diagram showing the proper motions of 
our spectroscopically confirmed T dwarfs with good proper mo- 
tion estimates (see section [3}. Shown are the different spectral 
types T4-6, T7-8 and T8+ as squares, triangles and circles, re- 
spectively. ULAS J1459+0857 and Wolf 940B are also highlighted. 



We searched around each of the 19 T dwarfs out to an 
angular separation corresponding to 20,000 AU, at the es- 
timated minimum distance of each T dwarf. We choose a 
separation limit of 20,000AU in order to be sensitive to 
the detection of both widely separated main-sequence (MS) 
and white dwarf (WD) companions. It is fairly common 
to find BD+ main-sequence star binaries with sepa rations 
of ~5000AU jGizis et al.l l200ll ; iPinfield et aHl2006ft . How- 
ever a white dwarf companion could have even wider sep- 
arations, when one considers any outward migration that 
would have occurred during the post-main-sequence mass- 
loss phase. We suggest that the outward migration would 
likely be up to a factor of ~4, as the initial and final sepa- 
ration of a low mass binary companion is directly rela ted to 
the change in mass of the host st ar i.e. Mms / Mwd (|Jeansl 
ll924l ; IZuckerman fc Becklinlll987| ). To illustrate this we con- 
sider a white dwarf of mass ~0.65 Mq (roughly the mean of 
the white dwarf mass distribution). The progenitor mass 
would be ~2.7Mq (from the initial-final mass relations o f 
iDobbie et al]|2006l ; ICatalan et aLll2008l ; iKalirai et aTll2008r i , 
such that Mms /Mot ~ 4. Thus, for BD + main-sequence 
binaries separated by ~5000AU, we could expect their final 
separation to be up to ~ 20,000AU. 

For each common proper motion companion candidate, 
we assumed a distance that was the same as the esti- 
mated distance of the relevant T dwarf, and thus placed 
companion candidates on an absolute B magnitude (Ms) 
vs B — R colour-magnitude diagram. Their positions were 
compared with those expected for main-sequence stars 
and white dwarfs, follow ing the approaches des cribed in 
IClarke et all |2009) and iDav- Jones et all [|2008n . respec- 
tively. We then compared their SDSS colours with respect 
to stell ar populations, includ ing main-se quence stars (Hip - 
parcos; IPerrvman et atll 19971). M dwar fs (|West et alj|2004 ). 
K subdwarfs jYong fc Lambertl 120031') and w hite dwarfs 
!|McCook fc Sionl 1 19991 . lEisenstein et al l I2006T ). We found 
that only one main-sequence companion candidate was 
identified, Wolf940 which had previously been identified 
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serendipitously by iBurningham et al l (|2009h and will not 
be discussed further. We also identified five candidate white 
dwarf companions to the T dwarfs in our sample, that are 
common proper motion (to within the uncertainties) and 
are consistent with the white dwarf sequence (see Fig. [3] and 
[4| if assumed to be at the same distance as their T dwarf 
companions. 

The brightest of our white dwarf candi dates appears in 
the L SPM catalogue (LSPM J1459+0851, iLepine fc Sharal 
I2005D as a high proper motion object, although it has 
not been previously studied spectroscopically. The T dwarf 
associated with this object is UL AS J 1459+0857, which 
has been spectroscopica lly typed as a T4.5±0.5 dwarf 
l|Burningham et al.l l2010). This pair is highlighted in Fig. [3] 
andU 

3.2 Possible contamination 

To assess the possibility that our candidate white dwarf com- 
panions may be dominated by high velocity background ob- 
jects, such as metal-poor, halo K-subdwarfs (which could 
populate the same colour and proper motion space that we 
search in this work), we have derived space motion esti- 
mates assuming tha t our candid a tes ar e K subdwarfs. We 
used the relations of llvezic et al.l (2008) to estimate an ab- 
solute r' magnitude (M r /) from g'-i' colour, for a metallic- 
ity range (for sub-dwarfs) of —0.5 and —2.5, and thus ob- 
tained distance constraints applicable if these objects are 
sub-dwarfs. Assuming for simplicity that they have a zero 
radial velocity, we then estimated space motions to assess 
potential halo membership. All except one candidate would 
have space motions of 2,000-12,000 kms -1 , which are thus 
not consistent with a galactic halo population. The excep- 
tion is close enough such that its kinematics are consistent 
with a background sub-dwarf member. Indeed this is one of 
the widest separated candidates (400 arcsec), and the level 
of background contamination for such separation (and vol- 
ume) app roaches 1, ev en for the low density halo luminosity 
function l|Gouldll2003h . However, despite one possible halo 
contaminant, the majority of the candidates cannot be ex- 
plained by such contamination. 

A more likely source of error potentially leading to 
misidentification amongst the candidates is proper motion 
uncertainty, since the ratio of proper motion to proper mo- 
tion uncertainty for some of the candidates is in the 3-5 
range, and one would thus expect that some fraction of the 
sample have proper motions that have scattered to larger 
values. Given this likely source of contamination we choose 
not to present details of all five of our companion candi- 
dates at this stage. We prefer to first establish their nature 
through spectroscopic study, and in that way confirm them 
(or otherwise) as genuine white dwarf companions. We have 
so far only obtained good spectroscopy for both components 
of the binary containing J1459+0851 and thus focus the re- 
mainder of this paper on this one system. 

4 SPECTROSCOPIC OBSERVATIONS OF 
LSPM J1459 + 0851 

Spectroscopic observations of LSPM J1459 + 0851 were ob- 
tained with FORS2 on the Very Large Telescope on 2009 



May 15 and 21, with Directors Discretionary Time in pro- 
gram 282.C-5069(A). We used the long slit mode in the op- 
tical wavelength range 3300-8000A with a dispersion of 50 
and 55A/mm respectively, for the ranges 3300-6210A (cor- 
responding to the B grism) and 5120-8450A (corresponding 
to the RI grism), giving a resolution of R~1200. Three inte- 
grations of 600s were taken, giving a total exposure time of 
30 minutes for the B grism and two integrations of 360s, to- 
taling 12 minutes in the RI grism. Sky flats, arc frames and 
the spectra of a DC white dwarf as well as a standard F- 
type star and DA white dwarf were taken during the same 
night at a similar airmass to the target so as to provide 
wavelength, flux and telluric calibrations. 

Standard iraf routines were used to reduce the spectra 
including flat fielding and cosmic ray removal. The spectra 
were extracted with APALL, using a chebyshev function to fit 
the background and a third order legendre function to trace 
the fit to the spectrum. The wavelength calibration was done 
using the spectrum from HgCdHeAr and HgCdHeNeAr arc 
lamps for the B and RI grisms respectively, and using iden- 
tify to reference the arc lines, along with the DISPCOR rou- 
tine to correct the dispersion of the spectrum. The resulting 
spectra of both LSPM J1459 + 0851 and the standard were 
divided by the smooth spectrum of the DC white dwarf, 
which has no intrinsic spectral features, enabling correction 
for the instrumental response. The standard stars (one for 
each grism) were then used to flux calibrate the spectrum. 
The two spectra were then stitched together in the overlap- 
ping sections and normalised at 6000A. The final spectrum 
of LSPM J1459 + 0851 is shown in Fig. [S] 

It can be seen that some residual tellurics remain and 
are highlighted for reference. They do not however affect 
the subsequent analysis in any way, since they do not over- 
lap with features used directly to assess white dwarf prop- 
erties. The general spectral shape is quite blackbody-like, 
consistent with a w hite dwarf or p erhaps a very metal 
poor sub-dwarf (e.g. |jao et al"1l2008l ). However the overall 
strength of the Ha line and the peak of the blackbody- 
like continuum are only c onsistent with a relatively cool 
example of a white dwarf (IKilic et al.ll2006l ). We also com- 
pare LSPM J 1459 + 0851 to the spectra of three other very 
cool, hydrogen rich (DA) white dwarfs, W D0011-134, WD 
1330+ 015 and WD 0503-174, ta ken from iBergeron et~ai1 
(1992) and IBergeron et alj (fl9 93). They have correspond- 
ing T eff 's of 6000+150 K, 7450+200 K and 5230+140 K re- 
spectively, and are shown in Fig. [6] Although the spectra 
of LSPM J1459 + 0851 are noisier, it can be seen that the 
extent of the Ha feature is consistent with a cool, hydrogen 
rich white dwarf. We further assess its properties in more 
detail in the following sections. 

4.1 Synthetic Photometry & Fitting Procedure 

As previously noted, the spectrum shows a lack of strong 
H lines, which would be expected for a hotter more typical 
white dwarf. Photometric fitting of the full optical-NIR SED 
is thus optimal for determining the effective temperature of 
the white dwarf, using model fits to the SDSS+UKIDSS 
photometry and assuming a distance equal to that of the 
T dwarf companion. The spectra are consistent with a cool 
white dwarf, showing no strong features in the spectrum 
blue-ward of 6000A and just a hint of Ha. We performed 
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Figure 3. A colour-magnitude diagram of white dwarfs from McCook & Sion (1999) with known parallax (crosses). Photometry is on 
the SuperCOSMOS system. Overplotted are model cooling tracks (see main text) for white dwarf masses of 0.5, 0.7 and 1.2 Mq (dotted, 
dashed and dot-dashed lines respectively). Also overplotted is our white dwarf selection region (two solid lines), along with our candidate 
white dwarf companions (large stars). LSPM J1459+0851 is circled for reference. 
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Figure 4. A two colour diagram in SDSS colours showing populations of main-sequence stars (blue points), M dwarfs (orange trian- 
gles), white dwarfs (green asterisks) and K subgiants (red squares), with our white dwarf candidates overplotted as black stars. LSPM 
J1459+0851 is circled for reference. 
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Figure 5. The optical spectra of LSPM J1459 + 0851, normalised 
at 6000A. 
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Figure 6. Optical spectrum in the region of Ha for the white 
dwarf LSPM J1459 + 0851. For comparison the spectra of three 
similar cool, hydrogen rich, magnetic white dwarfs, WD 0011- 
134 (B~16.7MG), WD 1330+015 (B~7.4MG) and WD 0503-174 
(B~7.3MG) are also shown, with magnetic field strength decreas- 
ing from top to bottom. 



a fit to the available photometry of the white dwarf (from 
SDSS, USNO and UKIDSS) using the atmospheric model 
codes of Berge r on et al., which are describ ed at length in 
Bergeron et al.l (1 19951 . with updates given in Ber geron et al.l 



200 li 120051 )." These models assume local thermodynamic 



equilibrium and allow energy transport by convection and 
can be calculated with arbitrary amounts of hydrogen and 
helium. Sy nthetic color were obtained using the procedure 
outlined in iHolberg fc Bergeron! (l2006h bas ed on the Vega 
fluxes taken from lBohlin fc Gillilandl l|20oi ). 

The method us ed to fit the photomet ric data is similar 
to that described in lBergeron et all (|200ll ). which we briefly 
summarize here. We first transform the magnitudes in each 
bandpass into observed average fluxes /™ using the following 
equation 



-2.5 log / A m + c„ 



where 



fV 

,/A 



J7 fxSm(X)d\ 



(1) 



(2) 



The transmission functions 5* m (A) alo ng with the constants 
Cm for each bandpass are described in lHolberg fc Bergeronl 
(2006) and references therein. To make use of all the photo- 
metric measurements simultaneously, we convert the magni- 
tudes into observed fluxes using equation fl} and compare 
the resulting energy distributions with those predicted from 
our model atmosphere calculations. Thus, we obtain a set 
of average fluxes /™, which can now be compared with the 
model fluxes. These model fluxes are also averaged over the 
filter bandpasses by substituting fx in equation ((2} for the 
monochromatic Eddington flux H\. The average observed 
fluxes /™ and model fluxes H x n , which depend on T e g , \ogg 
and 7V(He)/iV(H), are related by the equation 



fT = 4n(R/Dfm 



(3) 



where R/D is the ratio of the radius of the star to its dis- 
tance from Earth. Our fitting procedure relies on the nonlin- 
ear least-squares method of Levenberg-Marquardt, which is 
based on a steepest descent method. The value of \ 2 is taken 
as the sum over all bandpasses of the difference between both 
sides of equation (|3]), properly weighted by the correspond- 
ing observational uncertainties. In our fitting procedure, we 
consider on ly T e g and the s olid a ngle free parameters. As 
discussed in Bergeron et al.l (|200ll ). the energy distributions 
are not sensitive enough to surface gravity to constrain the 
value of log g, and thus for white dwarfs with no parallax 
measurement, as is the case here, we simply assume loggr 
= 8.0, which is consistent with the distance estimate of the 
T dwarf companion. Our best fit for a pure-H atmosphere 
arises from a T e s = 5535±45 K, and is shown in Fig. [7] 

4.2 A magnetic spectrum 

The shape and weakness of the Ha line gives a poor fit to 
a basic 5535 K model spectra, and we thus investigated the 
possibility that the white dwarf could be magnetic with the 
Ha line affected by Zeeman splitting. The line opacity was 



1 The synthetic colors can be obtained at 



http://www.astro.umontreal.ca/~bergeron/CoolingModels/ 
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calculated as the sum of the individual resonance-broadened 
Zeeman components. The line displacements and strengths 
of the Zeeman c omponents of Ha are taken from the tables 
of iKemicl (jl974), and the total line opacity is normalised to 
that resulting from the zero-field solution. The specific inten- 
sities at the surface, I(v,/j,,t v = 0), are obtained by solving 
the radiative transfer equation for various field strengths 
and values of [i (/x = cos#, where 6 is the angle between 
the angle of propagation of light and the normal to the sur- 
face of the star). In doing so, the polarization of the radi- 
ation is neglected as we are mainly interested in the total 
monochromatic intensity. The effect of the magnetic field on 
the continuum opacity is also neglected. The emergent spec- 
trum is then obtained from an integration over the surface 
of the star {H v oc J I^fidfi) for a particular geometry of the 
magnetic field distribution. We note that in this procedure, 
limb darkening is explicitly taken into account because of 
the integration over \i. 

We use the offset dipole model to model the magneti c 
field of the white dwarf l|Achilleos fe Wickramasingheil 19891 ) . 
In this model, the magnetic field is generated by a dipole. 
At the surface of the star (of radius unity), the strength of 
the magnetic field is simply given by 



B = 0.5B d (3cos 2 6>+ 1) 



1/2 



(4) 



where Bd is the dipole field strength, and 6 is the standard 
angle in spherical coordinates (9 = at the pole) . For a given 
value of Bd, the flux received at the Earth will also depend 
on the viewing angle i between the dipole axis and the line of 
sight (i — for a pole-on view). However, in this particular 
model, the dipole is also offset from the center of the star 
in an arbitrary direction. To simplify the calculation, we 
assume the dipole is offset parallel to the dipole axis. In this 
case, the value of the offset is measured from the center of 
the star and is denoted a z (in units of stellar radius). Note 
that with the offset dipole models, the value of the dipole 
field strength, Bd, is no longer equal to the value of the polar 
field strength. 

We computed a series of synthetic spectra based on a 
model atmosphere of pure-H with T c g — 5535 K and log g 
= 8.0 while varying Bd, i, and a z . We display in Fig. [7] 
the model which best reproduced the observed line profile 
with B d = 2.0 MG, i = 45° and a z = -0.20. It should 
be noted that varying the inclination angle i produced 
only slight var i ations in the line profiles and as stated in 
iBergeron etafl <| 19921 ). it is not possible to constrain i from 
observed line profiles alone. The zeeman-split Ha line in 
LSPM J1459 + 0851 can be compared to the same feature 
in other cool white dwarfs of higher magnetic field strengths 
in Fig [5] 



4.3 White dwarf age and progenitor lifetime 

We start by considering the simple case of a non- 
magnetic white dwarf for a T e g and log g of 5535 K 
and 8.0, for which we calculate a mass of 0.585 Mq 
l|Fontaine. Brassard fc Bergeronl 1200 ll ). The corresponding 
cooling age of a 0.585 Mq white dwarf, with a T e s = 
55 35 K was then calculated as 3 Gyr us ing the isochrones 
of iFontaine. Brassard fe "B ergeron (200lj). The total age of 
the white dwarf comprises of both the cooling time and its 
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Figure 7. Top: Fit of the energy distribution with pure-H mod- 
els. The observed ugriz and JHK fluxes are represented by error 
bars while the model fluxes averaged over the filter bandpasses 
are indicated by filled circles. The model monochromatic fluxes 
are shown by a solid line. Bottom: Normalised spectrum near 
Ha with the synthetic profiles interpolated at the parameters ob- 
tained from the energy distribution fits, assuming a pure-H at- 
mospheric composition with a model fit at T e ff = 5535 ±45K and 
logg = 8.0 with an offset dipole magnetic field computed with 
the parameters indicated in the figure. 



progenitor lifetime on the main-sequence. In order to con- 
strain the main-seque nce lifetime we acces s ed the Initial- 
final m ass relations of iFerrario et ail (|2005l ) , ICatalan et al.l 
l|2008h and lKalirai et al l (|2008l )~ to estimate a likely, initia- 
mass constraint for the main-sequence p rogenitor star of 
1.50-1 .75 Mq. We then used the tracks of lLachaume et aU 
(1999) to estimate the main-sequence lifetime for stars of 
such mass as 1.8-3.0 Gyrs. It should be noted however that 
the model tracks converge for masses <2Mq, and as a re- 
sult the ages of these objects can be largely uncertain (up 
to 10 Gyr). 

The strong magnetic field present also provides an ad- 
ditional factor to consider when assessing age. The origin of 
such strong white dwarf magnetic fields is not fully under- 
stood but is thought to have arisen in one of two favored 
scenarios: 

1. From a single star. The magnetic field is thought 
to derive from a massive, magnetic progenitor of ~1.5- 
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8Mq, typically an Ap or Bp star. The magnetic field 
is then maintained through the main-sequence evolu- 
tion to the white dwarf phas e by flux conservation 
|Wickramasinghe fc Ferrarioll200ol) . 

2. From the merger of two stellar cores in a common en- 
velope event or the merger of two degenerate objects. Dur- 
ing the common envelope (CE) phase the orbits of the two 
cores spiral in closer together through frictional forces caus- 
ing differential rotation, which coupled wit h convection in 
the c ores, creates a stellar-magnetic dynamo l|Tout fc Pingld 
1992). Close, but separated cores form CVs but some cores 
coalesce and cool to form a magnetic white dwarf. It may 
also be possible that two very close white dwarfs emerg e 
from the CE phase, such as G62-46 (|Bergeron et al.lll993l ). 
where one component is highly magnetic. 

In general it is observed that magnetic white dwarfs 
have larger masses than the mo re typical non-magnetic 
white dwarfs (e.g. iLiebertl [l988). There are two possi- 
ble hypothesis to explain this. Firstly, in accordance with 
the favoured scena rio for the formation of isolated mag- 
netic white dwarfs l|Wickramasinghe fc Ferrarid .2005'), the 
progenitor was more massive, leading to a massive white 
dwarf. In this case the magnetic field has no effect during 
the progenitor evolution. Secondly the effect of the mag- 
netic field has an impact on the stellar evolution, such 
that it could inhibit mass loss l|Wickramasinghe fc Ferrarid 
2000), leading to a more massive core and a longer progen- 
itor lifetime. If we consider the possibility that our white 
dwarf could be of higher mass, for example 0.8M q (the 
mean mass of a highly magnetic white dwarf; iKawka et al.l 
120071 ) t hen the cooling time would be lo nger, around 
6 Gyrs l|Fontaine. Brassard fc Bergeron! l200ll ) . In this case 
the progenitor (if a single star) would be around 3.5 Mq 
l|Catalan et al.ll2008f ) and the progenitor main-sequence life- 
time would be 0.3 Gyr (assuming normal models). However 
there is also evidence that many magnetic white dwarfs have 
masses closer to the pea k of the non-magnetic mass distri- 
bution (|Tout et al.ll2008l ). 

As we cannot know which scenario is responsible for 
the observed magnetic field, nor can we measure more ac- 
curately the mass of the white dwarf, we do not know the 
effects this may have had on the main-sequence evolution. In 
any case both scenarios for a magnetic and a non-magnetic 
white dwarf results in ages greater than 4.8 Gyrs (the cooling 
age plus the minimum progenitor lifetime of a non-magnetic 
white dwarf) and we thus choose to adopt this as the mini- 
mum total age for LSPM J1459 + 0851. 



Table 1. Parameters of the white dwarf LSPM J1459 + 0851. 



Parameter 



Value 



LSPM J1459+0851 - 
BOUND SYSTEM? 



ULAS J1459+0857 A 



In order to determine if this new system is a bonafide bi- 
nary system, we have statistically assessed the likelihood 
that two such objects could be a line-of-sight association 
with photometry and proper motion consistent with bina- 
rity by random chance. We first calculated the total re- 
gion of sky around our T dwarf corresponding to the cover- 
age encompassed by the projected line-of-sight separation 
of the white dwarf. We then combined this with the T 
dwarf distance constraint (43-69pc), allowing for the pos- 
sibility that the T dwarf might itself be an unresolved bi- 



RA 
DEC 
SDSS 'u 
SDSS 'g 
SDSS 'r 
SDSS 'i 
SDSS 'z 

SuperCOSMOS B 
SuperCOSMOS R 
SuperCOSMOS / 
USNO B 
USNO R 
USNO I 
UKIDSS Y 
UKIDSS J 
UKIDSS H 
UKIDSS K 
fj, RA 
DEC 

T e B 

Mass 
WD age 



14 59 32.05 
+08 51 28.1 
20.74 ± 0.08 
19.50 ± 0.01 
18.99 ± 0.01 
18.76 ± 0.05 
18.71 ± 0.03 
~19.48 
~18.33 
~18.29 
~19.8 
~18.8 
~18.2 
18.14 ± 0.02 
17.90 ± 0.02 
17.65 ± 0.04 
17.79 ± 0.06 
-170 ±3 mas yr" 1 * 
-42 ±6 mas yr —1 * 
5535±45 K 

8.0 dex 
0.585 M 
>4.8Gyr 



* USNO-B1 



nary (at a greater distance), to estimate a volume of sky 
in which white dwarfs might be line-of-sight contamina- 
tion. We then used the number density of white dwarfs 
(e.g. ISchroder. Pauli fc Napiwotzkil [2004 ) to estimate that 
we would expect only 0.00368 white dwarfs in this volume 
of space. 

To factor in the probability that two objects might 
have a common proper motion at the level of our measure- 
ments, we downloaded a magnitude-limited sample (i? <21, 
the same as our initial selection; see Section 13. 1|) of objects 
from the SuperCOSMOS Science Archive. We applied a limit 
to the proper motion uncertainty of <50mas/yr and required 
objects to lie in the colour range 1< B — R <3 (where we 
expect contaminant, main-sequence stars). This sample of 
10,360 sources was selected from within one degree of the 
T dwarf as to provide a representative sample of objects 
in the area of sky in which we find our binary system. We 
then placed these objects on a colour-magnitude diagram 
and selected only objects that occupied a region populated 
by main-sequence stars, when placed at the distance range 
estimated for the T dwarf. Of the 140 objects that were se- 
lected in this way, 13 had proper motion consistent (at the 
1 a level) with the T dwarf, leading to a probability of 0.092, 
that a contaminant star could have proper motion common 
with the T dwarf. 

The chance of finding a white dwarf with the 
same line-of-sight separation of LSPM J1459+851 from 
ULAS J1459+0857, where both components share a com- 
mon proper motion and are consistent with lying at the 
same distance as that estimated for the T dwarf is thus 
0.00368 x 0.092 = 0.0003. We thus conclude that these ob- 
jects form a genuine binary system. Since we searched for 
companions to a total of 19 T dwarfs, we estimate the over- 
all chance of finding a spurious system in our sample is 
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Table 2. Parameters of the binary system. 



Parameter 


Value 


Separation on sky 


385 arcsec 


Estimated distance 


43 - 69 pc* 


Estimated line-of-sight 




separation 


16, 500 - 26, 500 AU* 


* Assuming the T dwarf is 


a singular or unresolved binary. 



0.0064, and that the systems identified are likely real bi- 
naries. Properties of the new white dwarf + T dwarf sys- 
tem are given in Table. [21 and a finding chart is presented 
in Fig. [8] The wide separation of the system (16,500AU, 
assuming it is a singular object and not an unresolved bi- 
nary) is similar to the widest known BD+main sequenc e 
binary systems (e.g. iFahertv et alj|2009l . IZhang et al.ll2010h . 
Although prior to the post-main-sequence mass loss phase 
of the primary the separation would have been substantially 
less. Indeed we expect that the initial separation of the two 
components was a factor of ~4 closer (see Section ^) , in the 
region of~4100AU and more akin to the more common type 
of B D+main sequence binaries (<5000AU; IPinfield et all 
I2008T ). 



6 PROPERTIES OF ULAS J1459+0857 

In order to estimate the T cff of ULAS J 1459+0857 
we used spectral type-T^ff dete rminations from Ta- 
ble 6 of iGolimowski et all (|2004l ) as a guide. These 
determinations are for BDs with known parallaxes 



Tinnevl Il99fj; iLeggett et al.| 12002]; iDahn et al 
Tinncv, Burga sser fe Kirkpatrickl 120031 ; iKnapp et al 



2002 



2004 



Vrba et al.l |2004| ) and SEP constraints over a broad wave 



length range, for which reliable bolometric flux measure- 
ments and luminosities are thus available. The main source 
of uncertainty in these T c g values comes from a lack 
of strong age constraints and the resulting evolutionary 
model radii uncertain t ies (up to ~3 % for ages >0.1 Gyr s; 



Burrows et ail Il997l ; iBaraffe et ail Il998l ; IChabrier at al] 
2000). By considering the variety of T e s ranges calculated 



for the T4.5+0.5 dwarfs which have assumed a range of pos- 
sible ages from 0.1-10 Gyrs (± ~300K), we estimate that for 
an age range of 4-10 Gyrs ULAS J1459+0857 has T c g in the 
range 1200-1500K. 

We also used the Lyon Group COND models 
ijBaraffe et al.ll2~003h to estimate the physical properties of 
ULAS J 1459+0857, allowing for the possibility that it could 
be a single object or itself an unresolved binary, with a dis- 
tance in the range 43-69pc. We calculated absolute magni- 
tudes for ULAS J1459+0857 based on this distance range 
for ages 4-10 Gyrs. Then using a linear interpolation be- 
tween the model grid points we obtained mass and logp 
estimates in the range 0.06-0.072Mq and 5.4-5.5 dex respec- 
tively, assuming solar metallicity. We also consider the evo- 
lutionary models of Burrows et al. dBurrows et al.l Il997 



IBurrows et~aH l200ll ; iBurrows. Sudarskv fc Hubenvl l20oi 
to estimate a mass of 0. 064-0. 075Mq if the T dwarf is 
actually metal poor ([M/H] -0.5dex), which is a simi- 
lar value to that of the solar metallicity COND models. 



Both models also indicate a high gravity (logy = 5.5) for 
the observable J — K colours and temperature range of 
ULAS J1459+0857. We also compare the optical + NIR 
colours of ULAS J 1459+0857 in comparison with other spec- 
troscopically confirmed T dwarfs from the UKIDSS LAS 
l|Burningham et "aT1l2010h . Fig. I 111 shows this compliment of 
T dwarfs on a series of colour-spectral type plots. Whilst the 
NIR colours of ULAS J1459+085 in general look fairly typ- 
ical for a T4.5 dwarf, there is some evidence of relative z— 
band enhancement, which could be the result of the old er, 
higher gravity nature of this object l|Pinneld et al.ll2008l ). 



6.1 Initial model testing with benchmark 
observations 

To provide a first-pass test of model at mosphere predictions , 
we used the dust-free COND models of lBaraffe et~ai] (|2003l ) 
to provide theoretically informed best-guess constraints of 
the physical properties of the T dwarf in our binary. We 
made comparisons between the observed T dwarf spectrum 
and synthesised spectroscopy for two values of T c h (1200 
and 1500 K), three values of logg (5.0, 5.25 and 5.5) and 
three metallicities (+0.3, 0.0 and —0.5, representing the 



range observed in the galactic disk; e.g. IValenti fc Fischer 



2005; Gray et al. 2006; Holm berg et al.ll2007l ; I Jenkins et al 



2008 



The model spectra were generated using the atmo- 
spheric rad iative transfer code, Phoenix (which is described 
in de tail in lHauschildt. Allard fc Baron 1 ll999l ; lAllard" et al.l 
l200lf ) which includes the upda ted water molecular opac- 
ity li st from Barber et al.l (120071) and new solar abundances 
from lAsphmd et al. ( 2009T ) . The model includes the effect 
of condensation in the chemical equilibrium but ignores the 
effects of dust opacities. Theoretical spectra were normalised 
to the observations in the peak of the J— band, and the re- 
sulting comparisons are shown in Fig. [9] and 1101 The long 
red dotted line, short dashed green and the long dashed blue 
lines represent the logy = 5.0, 5.25 and 5.5, respectively. 

Table. [3] gives a summary of the reduced chi-squared 
ixV) values calculated for these comparisons. These values 
range from ~1.75-3.5, and we note that none of them are 
close to 1.0. This is to be expected since there are known 
shortcomings in the models that introduce differences signif- 
icantly greater than the measurement uncertainties - e.g. in- 
complete methane opacities and a poor underst anding of the 
obser ved J— band brightnening around T3 (e.g. lKnapp et al] 
l2004h . As such, we only use these xt values as an additional 
quality indicator for the over-all model-observation compar- 
ison. 

Our analysis reveals that the best fitting atmospheric 
model has an effective temperature of 1200K, a sub-solar 
metallicity and a high gravity of logg =5.5. However, com- 
parison by eye shows clearly that this fit across the whole 
spectral range is not particularly good, especially around 
the peak flux regions of the H— and K— bands. This best 
Xt value presumably comes from the better fit to the J— 
band, compared to the other models, however a better fit 
(by eye) to the overall profile of the spectra comes from a 
1500K, solar metallicity, high gravity model. This suggests 
that ULAS J1459+0857 probably has solar to slightly sub- 
solar metallicity and high gravity. This is instructive at least 
that the observed properties of ULAS J1459+0857 appear to 
be in general agreement with those predicted by evolution- 
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Figure 8. A UKIDSS J band finder chart showing the position of the T dwarf (square) and the white dwarf (circle), the scale of the 
image is 8' x 3.5'. 



Table 3. xt n * s to model spectra. 



Table 4. Parameters of ULAS J1459 + 0857. 



log 9 






[M/H] 








+0.3 


0.0 


-0.5 


T e g= 


1200 K 








5.00 




1.94 


2.11 


2.38 


5.25 




2.06 


2.58 


1.82 


5.50 




1.97 


2.07 


1.76 


T e s = 


1500 K 








5.00 




3.50 


2.47 


2.55 


5.25 




2.89 


2.03 


2.62 


5.50 




3.09 


1.96 


2.88 



ary models. While the models clearly have some shortcom- 
ings they appear to be making progress, such that model fit 
properties for a mid T dwarf of age > 4Gyrs seem to be 
broadly consistent with the benchmark fit properties. 



7 SUMMARY 

This is the first discovery of a T dwarf + white dwarf binary 
system, and an example of an evolved, high gravity brown 
dwarf. During the main sequence phase of the primary its 
separation would have been similar to the bulk population of 
BD+main sequence binaries, although the separation must 
have grown significantly (to its current value) during the 
post-main-sequnce mass loss phase of the primary. 

The white dwarf provides vital age constraints for the 
binary system from it's cooling age combined with a mini- 
mum estimate of it's main sequence progenitor lifetime, and 
the resulting minimum age of 4.8 Gyrs for the T dwarf allows 
a robust constraint on its surface gravity of log g = 5.4 — 5.5. 
As such ULAS J1459+0857 is a representative old, high 
gravity benchmark BD. Comparison with the bulk popual- 
tion of UKIDSS T dwarfs shows some indication that z— 
band flux enhancement may be an observational character- 
istic of older high gravity mid T dwarfs, which could be a 
useful factor in attempts to understand the formation his- 



Parameter Value 

HA 14 59 35.25 

DEC +08 57 51.20 

Distance 43 - 69pc* 1 

SDSS z' 21.17+0.15 

UKIDSS Y 19.24+0.06 

UKIDSS J 17.93+0.02 

UKIDSS H 17.94+0.05 

UKIDSS K 17.92+0.08 

UFTI J 17.98 ± 0.04 

UFTI H 17.93+0.04 

UFTI K 18.04+0.03 

|U RA -149 ±33 mas yr -1 

H DEC -45 ±33 mas yr" 1 

Spectral Type T4.5 ±0.5 

Mass 0.060 - 0.072 M * 2 

Mass 0.064 - 0.075 M * 3 

T cff 1200 - 1500 K 

logp 5.4 — 5.5 dex 

41 Assuming the T dwarf is a singular or unresolved binary. 
* 2 From the Lyon group COND evolutionary models. 
* 3 From the Burrows evolutionary models. 



tory of sub-stellar objects. And more generally, this T dwarf 
can contribute to our understanding of substellar properties 
by providing a useful test-bed for theoretical model atmo- 
spheres. 

This system is an example of how wide BD binary 
companions to white dwarfs make good benchmark objects, 
which will help test model atmospheres, and may provide 
indep endent means to ca librate BD properties of field ob- 
jects (|Pinfield e t al. 2006). This is the first of our candidate 
systems and we expect to find many more as we mine more 
UKIDSS sky, and with our ongoing efforts to search still 
deeper by combining UKIDSS BDs with SDSS for white 
dwarfs. 

Further observations that would have clear benefits for 
this system include a parallax measurement of either compo- 
nent to yield an accurate distance, and fuller spectral cover- 
age of the T dwarf to constrain its mid-infrared and optical 
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Figure 10. Spectral model comparisons to ULAS J1459+0857 (black solid line) for T off =1500K and [M/H]= +0.3, 0.0 and -0.5, from 
top to bottom, with logg = 5.0, 5.25 and 5.5 as a long red dotted line, short dashed green and the long dashed blue lines, respectively. 
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Figure 11. Colour against T dwarf spectral type for UKIDSS T dwarfs jBurningham et al]|2010h . ULAS J1459+0857 is shown as a 
large, filled red circle. 



spectral morphology. In general, and particularly in the op- 
tical z— band, it would be desirable to assess in more detail 
spectroscopic features and trends that may be sensitive to 
higher surface gravity and age. A parallax distance com- 
bined with accurate knowledge of the T dwarf bolometric 
flux would offer significant improvements on the T e g con- 
straints for this object fe.g. iBurningham et aLlfeoOSlI ). 

An accurate distance would also facilitate greatly im- 
proved radius and mass constraints for the white dwarf, and 
thus a better constraint on the cooling age and progenitor 
lifetime. More detailed studies of both binary components 
are clearly crucial to maximize the effectiveness of the bench- 
mark BD component. 
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